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Proposed silicon-based quantum-computer architectures have attracted attention because of their 
promise for scalabihty and their potential for synergetically utilizing the available resources associ- 
ated with the existing Si technology infrastructure. Electronic and nuclear spins of shallow donors 
(e.g. phosphorus) in Si are ideal candidates for qubits in such proposals because of their long spin 
coherence times due to their limited interactions with their environments. For these spin qubits, 
shallow donor exchange gates are frequently invoked to perform two-qubit operations. We discuss 
in this review a particularly important spin decoherence channel, and bandstructure effects on the 
exchange gate control. Specifically, we review our work on donor electron spin spectral diffusion 
due to background nuclear spin flip-flops, and how isotopic purification of silicon can significantly 
enhance the electron spin dephasing time. We then review our calculation of donor electron ex- 
change coupling in the presence of degenerate silicon conduction band valleys. We show that valley 
interference leads to orders of magnitude variations in electron exchange coupling when donor con- 
figurations are changed on an atomic scale. These studies illustrate the substantial potential that 
donor electron/nuclear spins in silicon have as candidates for qubits and simultaneously the consid- 
erable challenges they pose. In particular, our work on spin decoherence through spectral diffusion 
points to the possible importance of isotopic purification in the fabrication of scalable solid state 
quantum computer architectures. We also provide a critical comparison between the two main pro- 
posed spin-based solid state quantum computer architectures, namely, shallow donor bound states 
in Si and localized quantum dot states in GaAs. 

PACS numbers: 03.67.Lx, 71.55.Cn 
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I. INTRODUCTION 



During the past decade, the study of quantum computing and quantum information processing has generated 
widespread interest among physicists from areas ranging from atomic physics, optics, to various branches of condensed 
matter physics ■ The key thrust behind the rush toward a working quantum computer (QC) is the development of 
a quantum algorithm that can factorize large numbers exponentially faster than any available classical algorithm j^]- 
This exponential speedup is due to the intrinsic quantum parallelism in the superposition principle and the unitary 
evolution of quantum mechanics. It implies that a computer made up of entirely quantum mechanical parts (qubits), 
whose evolution is governed by quantum mechanics, would be able to carry out prime factorization of large numbers 
that is prohibitively time-consuming in classical computation, thus revolutionizing cryptography and information 
theory. Since the invention of Shor's factoring algorithm, it has also been shown that error correction can be done to 
a quantum system so that a practical QC does not have to be forever perfect to be useful, as long as quantum 
error corrections can be carried out on the QC. These two key mathematical developments have led to the creation 
of the new interdisciplinary field of quantum computation and quantum information. 

Many physical systems have been proposed as candidates for qubits in a QC. Among the more prominent exarnples 
are electron or nuclear spins in semiconductors |^ Q , includin g el ectron spin in semiconductor quantum dots 0, ll| 
and donor electron or nuclear spins in semiconductors ITollTH. The donor-based QC schemes are particularly 
interesting because all donor electron wavefunctions in a semiconductor are identical, and because doping makes a 
natural connection between quantum mechanical devices and the more traditional microelectronic devices: Doping 
in semiconductors has had significant technological impact for the past fifty years and is the basis of the existing 
microelectronics technology. As transistors and integrated circuits decrease in size, the physical properties of the 
devices are becoming sensitive to the actual configuration of impurities 0|. In this context, the important proposal 
of donor-based silicon quantum computer (QC) by Kane @, in which the nuclear spins of the monovalent "^^P impurities 
in Si are the qubits, has naturally created considerable interest in revisiting all aspects of the donor impurity problem 
in silicon, particularly in the Sii'^^P system. 

In principle, both electron spin and orbital degrees of freedom can be used as qubits in semiconductor nanostructures. 
For example, electron orbital dynamics is quantized into discrete "atomic-like" levels (on meV energy scale with a 
Bohr radius of the order of 10 nm) in semiconductor quantum dots, and two such quantized quantum dot levels could 
form the quantum two-level system needed for a qubit JT^ . A great advantage of such orbital (or equivalently, charge) 
qubits is that qubit-specific measurements are relatively simple since one is essentially measuring single charge states, 
which is a well-developed experimental technique through single-electron transistors (SET) or equivalent devices [l3| . 
A major disadvantage of solid state charge qubits is that these orbital states are highly susceptible to interactions 
with the environment (which, in particular, contains all the stray or unintended charges inevitably present in the 
device), and the decoherence time is generally far too short (typically picoseconds to nanoseconds) for quantum error 
correction to be useful. A related problem is that inter-qubit coupling, which is necessary for the implementation 
of two-qubit gate operations essential for quantum computation, is often the long-range dipolar coupling for charge 
qubits. This makes it difficult to scale up the architecture, since decoherence grows with the scaling-up as more and 
more qubits couple to each other via the long-range dipolar coupling. Additional decoupling techniques have to be 
applied to ensure selective gate operations and to control decoherence due to long range couphngs |l^. Since 
scalability is thought to be the main advantage of solid state QC architectures, little serious attention has so far been 
paid to orbital qubit based quantum computation in semiconductor nanostructures due to its probably unscalable 
decoherence and entanglement properties. 

Spin qubits in semiconductor nanostructures have complementary advantages (and disadvantages) compared with 
charge qubits based on quantized orbital states. A real disadvantage of spin qubits is that a single electron spin is 
difficult to measure, although there is no fundamental principle against the measurement of a Bohr magneton. The 
great advantage of spin qubits is the very long spin coherence times, which even for electron spins can be milliseconds 
(microseconds) in silicon (GaAs) at low temperatures. This six orders of magnitude coherence advantage in spin 
qubits over charge qubits has led to electron spin qubits in GaAs quantum dots and in P donor levels in silicon (as 
well as SiCe quantum structures) as the QC architectures of choice for the solid state community. In addition to 
the coherence advantage, spin qubits also have a considerable advantage that the exchange gate, which provides the 
inter-qubit coupling, is exponentially short-ranged and nearest-neighbor in nature, thus allowing precise control and 
manipulation of two-qubit gates. There is no fundamental problem arising from the scaling-up of the QC architecture 
since exchange interaction couples only two nearest-neighbor spin qubits independent of the number of qubits. In this 
review we provide a brief perspective on spin qubits in silicon with electron spins in shallow P donor levels in Si being 
used as qubits. 

Although experimental progress in semiconductor-based solid state QC schemes has been slow to come during 
the past five years, they are still often considered promising in the long term because of their perceived scalability 
advantages. After all, the present computer technology is based on semiconductor integrated circuits with ever smaller 
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feature size. Therefore, semiconductor nanostructure-based QC architectures should in principle be scalable using 
the existing microelectronics technology. However, it still remains to be demonstrated whether and how the available 
(classical) semiconductor technology can help the scaling up of a quantum coherent QC architecture. For the spin 
qubits in silicon, for example, the key issues include clarifications of spin quantum coherence properties in the solid 
state environment, physical approaches to manipulate and entangle spins, fabrication of devices with atomic-scale 
precision, and measurement of single spins. In the following, we review our work on two of these important issues: 
donor electron spin coherence and spin interaction in silicon. 



II. SPIN COHERENCE IN SILICON NANOSTRUCTURES 

Before the seminal concept of quantum error correction was introduced in 1995 Q, it was widely believed that 
quantum computation, even as a matter of principle, is quite impossible since all quantum states decohere due to 
interaction with the environment, and such decoherence was thought to be fatal to QC operations. Although the 
quantum error correction principle has shown that a certain degree of decoherence can be corrected in QC algorithms, 
one still has severe limits on the amount of tolerable decoherence for feasible QC operations. The question of 
decoherence is therefore of paramount importance for any form of quantum computing architecture, including the 
spin-based proposals. Estimates of exchange coupling and various analysis of adiabatic operation of the exchange 
gates in quantum dots IT] suggest that the exchange gate duration is unlikely to be shorter than 100 ps-adiabaticity 
is likely to be compromised for shorter pulse times leading to loss of qubit fidelity. Hence spin coherence times at 
least of the order of microseconds are necessary to satisfy the criteria of fault-tolerance (10'* reliable quantum gate 
operations during the coherence time Jj). Non-ideal situations encountered in real devices will probably require even 
longer coherence times. Here we argue that the possibility of using nuclear-spin-free samples of silicon through isotopic 
purification together with the small spin-orbit coupling in silicon makes Si one of the most promising host materials 
with respect to electron spin coherence. As we show below, both adiabaticity (i.e. "slow" gate operation) and fault 
tolerance (i.e. "fast" gate operation) are more easily satisfied in Si donor electron spins than in GaAs quantum dots. 

The simplest description of electron spin coherence is to consider characteristic time scales Ti and T2 [1^ , which are 
respectively defined as the decay times of spin magnetization parallel (the so-called longitudinal spin relaxation time 
Ti) and perpendicular (the so-called transverse spin dephasing time T2) to an external magnetic field. All processes 
contributing to Ti require energy exchange with the lattice (via phonon emission, for example). These Ti relaxation 
processes also contribute to the dephasing time T2, hence the inequality T2 < 2Ti 19]. However, often one finds 
T2 <C Ti due to the importance of dephasing mechanisms that do not involve energy relaxation, thus contributing 
exclusively to T2 decay. This is exactly the case for a Si:P donor electron spin. In this system Ti can be as long as 
10^ seconds |23] (measured for B = 0.3 T - we note that Ti oc B^^ at low temperatures due to phonon matrix element 
and phase space considerations ^^), while T2 detected by spin echo decay is of the order of 10~^ seconds The 
reason why Ti is so long in silicon nanostructures is two fold: (1) Spin flips mediated by spin-orbit coupling and the 
electron-phonon interaction are strongly suppressed for localized electrons |22j due to energy-momentum conservation 
constraints, and (2) spin-orbit coupling in silicon is significantly weaker than in other popular semiconductors 
(such as III-V compounds, where Ti is of the order of nanoseconds for conduction electrons Q). This observation, 
together with the fact that one should expect T2 = 2Ti in isotopically purified ^^Si (see below), indicates that silicon 
is one of the most attractive materials for coherent spin manipulation. 

It was discovered a long time ago :24!j that the dipolar fluctuation of lattice nuclear spins is responsible for electron 
spin phase fluctuations, leading to signiflcant suppression of T2, an effect usually denoted spectral diffusion. Here each 
nuclear spin produces a different hyperflne shift in the electron Zeeman frequency. Because the nuclei are coupled 
to each other via magnetic dipolar interaction leading to nuclear spin flip-flop transitions, the total hyperflne fleld 
at the localized electron fluctuates and leads to time-dependent noise in the electron spin Zeeman frequency and 
its consequent decoherence. Below we summarize a microscopic theory for this spectral diffusion effect which was 
proposed by two of us |2^ 12^ and recently verifled experimentally [2^ [281 ] . 

Silicon has three stable isotopes: ^^Si has nuclear spin quantum number / = 1/2 and a natural abundance equal to 
4.68%. ^^Si and ^°Si are spin-0 and hence do not contribute to spectral diffusion. Nuclear isotope engineering |27ll2q 
can systematically reduce the amount of ^^Si ("the isotopic purification"), thus raising an important theoretical 
question on how T2 depends on the fraction of ■^^Si present in the sample. (We note that T2 due to spectral diffusion 
becomes infinitely long as the fraction of ^^Si in the sample approaches zero.) To address this question, we consider 
a truncated Hamiltonian for the electron-nuclear spin evolution. 
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FIG. 1: Electron spin coherence time for a Si:P donor. 9 denotes the angle between the external magnetic field and the [100] 
crystallographic direction. The agreement between theory and experiment is within a factor of 3 (experimental data from 
Ref. IH). 
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Here Sz is the z-component of the electron spin operator, and 7s is the electron gyromagnetic ratio. I„ is the nuclear 
spin operator for a ^^Si isotope located at position Rn, An is the hyperfine shift produced by this nucleus, 



and bnm is the dipolar coupling between two such nuclei, 
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which has an important dependence with respect to the angle 0„,„ between the magnetic field and the vector linking 
the two nuclei R„ — R^. In Eq. ^ we neglect the off-diagonal hyperfine coupling {S+In- + h.c), which can be 
shown to play no role at magnetic fields higher than the Overhauser field {B > An ~ 10 — 1000 G in silicon [2^). 
The last term of Eq. is responsible for flip-flop processes between two nuclear spins, which can be studied using 
an inverse temperature expansion. The rate for a flip-flop process is then given by [26j 
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Here A„ 



\An — Am\/2 is the frequency shift felt by the electron when one flip-flop event takes place, while the 



linewidth for flip-flop is given by 
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Nuclear spin flip-flop can only take place if the nuclear system can rearrange itself to compensate for the energy cost 
A„„i- Hence pairs of nuclei satisfying A„m ^ Unm are inert. Since these are usually located close to the center of the 
donor (where An is larger), they form a frozen core which does not contribute appreciably to spectral diffusion. 

An approximate theory for spectral diffusion decay can be obtained by assuming that the most important processes 
contributing to nuclear dipolar fluctuation are uncorrelated flip-flop transitions with Markovian dynamics. The 
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FIG. 2: Dependence of Si:P electron spin coherence time (as measured by Hahn echo decay) with the fraction of '^^Si nuclear 
spins. This figure reveals isotopic purification as an effective tool for material optimization of T2. 



7r/2 — r — TT — T — echo envelope arising due to the stochastic fluctuation of one pair of nuclei n, m is then determined 
to be 

Vnm{2T) = i?-,?,C-2-r"-|r2^C0Sh(2Ti?„,„) 

+r„,„i?„m sinh (2ri?„,„) - A^„„|, (6) 

where Rnm — T^nm ~ ^nm- We then consider the contribution of all nuclear spins through an averaged product 

v{2r)= n [(l-/2) + /2^;„„(2r)] , (7) 

where / is the fraction of ^^Si nuclear spins present in the lattice. Equation Q allows realistic calculations of echo 
decay due to spectral diffusion without using any fitting parameters. All we need is an appropriate wave function 
for the electron (we use the Kohn-Luttinger state for a Si:P impurity) which in turn determines the hyperfinc shifts 
A„m and the corresponding fluctuation rates r„m (other constants such as 7s, 7/ are readily available from the 
NMR literature). For natural silicon, the calculated T2 is 2.5 times larger than the measured value j^J, a reasonable 
agreement in the spin relaxation literature @ . In addition we were able to predict an important orientation dependence 
for T2 (see Fig. ^ . As the direction of the external magnetic field is changed with respect to the silicon diamond 
structure To increases by as much as a factor of three. Recently, this prediction was verified experimentally for natural 
silicon [23, nil and for a ^^Si enriched sample |23| (Fig. 0). The calculated dependence of T2 on the fraction of ^^Si 
is shown in Fig. |21 establishing isotopic purification as an eflticient way to achieve the longest possible electron spin 
coherence times for QC use. Recently A. M. Tyryshkin and collaborators [285 performed echo spectroscopy in a sample 
with less than 0.005% ^^Si isotopes. The residual Ti > 60 ms was attributed to the lack of coherence in the applied 
ESR pulses, and not to the spectral diffusion process. 

Isotopic purification can also be used in germanium, where the active nuclear spin is ^^Ge (/ = 9/2, 7.73% natural 
abundance). On the other hand the important class of IITV semiconductors (GaAs, GaSb, InAs, InSb) has no known 
/ = nuclear isotopes. Because echo decay has never been measured for localize d sp ins in these materials, the role 
of theory becomes extremely important. Recently we predicted T2 ^ I — 100 /iS |2^|2^ for localized electron spins 
in III-V semiconductor quantum dots, where the development of single spin detection techniques promises a bright 
future for coherent spin manipulations pof . There are currently two proposals to reduce the effect of strong spectral 
diffusion associated with the absence of / = isotopes in III-V materials. The first relies on substantial nuclear 
polarization [s^l (very high nuclear polarization is required to suppress the flip-flop processes that lead to spectral 
diffusion p5l |. The qualitative dependence of T2 on nuclear polarization can be derived from the short r limit of 
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FIG. 3: Depicts the dependence of coherence time T2 with the nuclear spin polarization p. Very high nuclear polarization 
(p > 0.99) is required to enhance T2, a condition which is very hard to achieve experimentally. 



Eqs. © and |(7J. In this case we have (2^ 

L \' 

-\ (xp^pi(x{l-p^), (8) 

where the nuclear spin polarization is given hy p — p-\ — pi- This dependence implies T2 is enhanced significantly only 
if p is very close to 1 (Fig.O. For example, a nuclear polarization ofp = 99.5% increases T2 only by a factor of 5). The 
other approach to reduce spectral diffusion requires successive applications of tt pulses (Carr-Purcell-Meiboom-Gill 
sequence) at a rate proportional to the square of the nuclear spin quantum number l3l| . Both these proposals require 
considerable overhead, and make evident the advantages of silicon or germanium from the perspective of coherent 
spin manipulation. 

We make two brief final comments before concluding our discussion of spin decoherence. The first comment is that 
we believe the spectral diffusion mechanism to be the ultimate decoherence mechanism for localized electron spins in 
semiconductor nanostructures since it cannot be eliminated in a given sample by reducing the sample temperature 
as the nuclear spin energy scale is extremely small (sub-millikelvins). Thus, understanding spectral diffusion effects 
in spin-based solid state QC architectures is one of the most important physics issues in the context of electron spin 
dynamics in quantum computers. From this narrow perspective, silicon QC has significant advantages over GaAs QC 
due to the unique role of isotopic purification in silicon. The second comment is that our statement about the spin 
coherence time in GaAs nanostructures (or more generally, in III-V materials) being 'short' {T2 ~ 1 — 100 /is) should 
be understood in the proper context of a spectral diffusion comparison with the corresponding silicon coherence where 
T2 ^ 1 — 100 ms. In this context, recent claims in the literature of very long spin coherence times 1-100 ns) in 
GaAs for conduction electrons are in fact extraordinarily short for the purpose of spin quantum computation. 



III. DONOR ELECTRON EXCHANGE IN SILICON 



An important issue in the study of donor-based Si QC architecture is coherent manipulations of spin states as 
required for the quantum gate operations. In particular, two-qubit operations, which are required for a universal QC, 
involve precise control over electron-electron exchange 0, 0, llTl Is^ and electron-nucleus hyperfine interactions (for 
nuclear spin qubits) . Such control can presumably be achieved b y fa brication of donor arrays with accurate positioning 
and surface gates whose potential can be precisely controlled |33l Is^. Issl Is^ . However, electron exchange in bulk 
silicon has spatial oscillations [s^, HI] on the atomic scale due to valley interference arising from the particular six- fold 
degeneracy of the bulk Si conduction band. These exchange oscillations place heavy burdens on device fabrication 
and coherent control 38] , because of the very high accuracy and tolerance requirements for placing each donor inside 
the Si unit cell, and/or for controlling the external gate voltages. 

The potentially severe consequences of the exchange-oscillation problem for exchange-based Si QC architecture 
motivated us and other researchers to perform theoretical studies going beyond some of the simplifying approximations 
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in the formalism adopted in Ref. ^3§, and incorporating perturbation effects due to applied strainjs^ or gate fields 
[i^ l. These studies, all performed within the standard Heitler-London (HL) formalism essentially reconfirm 
the originally reported difficulties regarding the sensitivity of the electron exchange coupling to precise atomic-level 
donor positioning, indicating that they may not be completely overcome by applying strain or electric fields. The 
sensitivity of the calculated exchange coupling to donor relative position originates from interference between the 
plane-wave parts of the six degenerate Bloch states associated with the Si conduction-band minima. More recently 
[42] we have assessed the robustness of the HL approximation for the two-electron donor-pair states by relaxing the 
phase pinning at donor sites, which could in principle eliminate the oscillatory exchange behavior. Within this more 
general theoretical scheme, the floating-phase HL approach, our main conclusion is that, for all practical purposes, 
the previously adopted HL wavefunctions are robust, and the oscillatory behavior obtained in Refs . ISa . l39l liol persists 
(quantitatively) in the more sophisticated theory of Ref. ^3 ■ 

In what follows, we first review the main results leading to the exchange oscillation behavior qualitatively described 
above. We then consider two substitutional donors in bulk Si, and present a systematic statistical study of the 
correlation between the relative position distributions and the resulting exchange distributions. We also show that 
strain may partially alleviate the exchange oscillatory behavior, but it cannot entirely overcome it. 



A. Donor Electron Exchange in Relaxed Bulk Silicon 

We describe the single donor electron ground state using the effective mass theory. The bound donor electron 
Hamiltonian for an impurity at site Ro is written as 

Ti-a = Ti-sv + Ti-vo ■ (9) 

The first term, Tisv, is the single- valley Kohn-Luttinger Hamiltonian |43j |. which includes the single particle kinetic 
energy, the Si periodic potential, and the screened Coulomb perturbation potential produced by the impurity ion 

For shallow donors in Si, we use the static dielectric constant e — 12.1. The second term of Eq. lO, Hvo, represents 
the inter-valley scattering effects due to the presence of the impurity. 

The donor electron eigenfunctions are written in the basis of the six unperturbed Si band edge Bloch states 
(pfj^ = u^(r)e'''f' '' [recall that the conduction band of bulk Si has six degenerate minima (/i = 1, . . . , 6), located along 
the F— X axes of the Brillouin zone at |k^| ~ 0.85(27r/a) from the F point]: 

^^oir) = ^^F^(r-Ro)^^(r)e'''--("-^o). (11) 

The phases of the plane- wave part of all band edge Bloch states are naturally chosen to be pinned at Rq: In this 
way the charge density at the donor site [where the donor perturbation potential energy Eq. (|10|) is the smallest] is 
maximum, thus minimizing the energy for ipnoi'"^)- 

In Eq. I|ll|l . F^(r — Rq) are envelope functions centered at Rq, for which we adopt the anisotropic Kohn-Luttinger 
form, e.g., for fi = z, Fz{r) = exp{— [(x^ + y'^)/a'^ + /b"^] ^^^j/V -Ko^b. The effective Bohr radii a and b are variational 
parameters chosen to minimize Esv — (V'Rol^sv|'0Ro)i leading to a = 25 A, 6 = 14 A and Esv ~ —30 meV when 
recently measured effective mass values are used in the minimization [s^ ]. The periodic part of each Bloch function 
is pinned to the lattice, independent of the donor site. 

The TLsv ground state is six-fold degenerate due to the six-fold valley degeneracy of Si conduction band. This 
degeneracy is lifted by the valley-orbit interactions 44], which are included here in Tivoj leading to the nondegenerate 
(Ai-symmetry) ground state in (|ll|l . 

The HL approximation is a reliable scheme for the well-separated donor pair problem (interdonor distance much 
larger than the donor Bohr radii) Within HL, the lowest energy singlet and triplet wavefunctions for two electrons 
bound to a donor pair at sites R^ and Rs, are written as properly symmetrized combinations of ipn^ and i/'Rb [as 
defined in Ea. l(TT)l ] 

^tiri,r2) = , [V^R^(ri)VR^(r2) ± ^r Jri)^R^ (ra)] , (12) 

V2(l ± 5^) 
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FIG. 4: Exchange coupling between two phosphorus donors in Si along the indicated directions in the diamond structure. 
Values appropriate for impurities at substitutional sites are given by the squares. The dashed line in the _R||[100] frame is a 
guide to the eye, indicating that the oscillatory behavior may be ignored for donors positioned exactly along this axis. 

where S is the overlap integral and the upper (lower) sign corresponds to the singlet (triplet) state. The energy 
expectation values for these states, = (^'f give the exchange splitting through their difference, J = Et — Eg. 

We have previously derived the expression for the donor electron exchange splitting |33,E3j which we reproduce here: 

'^W^^E'^M-(R)cos(k^-k,).R, (13) 

where R = — Rb is the interdonor position vector and J7^^ (R) are kernels determined by the envelopes and are 
slowly varying [s^y^- Note that Eq. H13|l does not involve any oscillatory contribution from w^(r), the periodic 
part of the Bloch functions [iol li^ . The physical reason for that is clear from Hll|) : While the plane- wave phases of 
the Bloch functions are pinned to the donor sites, leading to the cosine factors in H13|l . the periodic functions are 
pinned to the lattice, regardless of the donor location. 

The exchange energy calculated from Eq. H13() for a pair of donors as a function of their relative position along the 
[100] and [110] crystal axis is given in Fig.0] This figure vividly illustrates both the anisotropic and the oscillatory 
iDchavior of J(R), which is well established from previous studies [s^, l38L l3^ liol ]. It is interesting to note that 
for substitutional donors with interdonor position vectors exactly aligned with the [100] crystal axis, the oscillatory 
behavior may be ignored in practice, as indicated by the dashed line in the figure. This behavior is qualitatively 
similar to the exchange versus donor separation dependence assumed in Kane's proposal 9], where the Herring and 
Flicker expression 45], originally derived for H atoms, was adapted for donors in Si. Therefore one might expect that 
reliable exchange gate operations would be possible if all donor pairs are exactly aligned along the [100] crystal axis. 

B. Nanofabrication aspects 

Aiming at the fabrication of a P donor array accurately positioned along the [100] axis, and taking into account the 
current state of the arts degree of control in substitutional P positioning in Si of a few nm [33, 34, 35, 36], we investigate 
the consequences of such interdonor positioning uncertainties {'^ a few nm) in the values of the corresponding pairwise 
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FIG. 5: Interdonor distance distributions for a target relative position of 20 lattice parameters along [100] (see arrows). The 
first donor is fixed and the second one "visits" all of the Si substitutional lattice sites within a sphere centered at the target 
position, with uncertainty radii (a)10 A, (b)20 A and (c)50 A. 



exchange coupling. We define the target interdonor position Rt along [100], with an arbitrarily chosen length of 
20 lattice constants (~ 108.6 A) indicated by the arrow in Fig. 01 The distributions for the interdonor distances 
R = [Ryi — Rb| when Ha is fixed and Rb "visits" all of the diamond lattice sites within a sphere centered at the 
target position are given in Fig. [51 Different frames give results for different uncertainty radii, and, as expected, 
increasing the uncertainty radius results in a broader distribution around the target distance. Note that the geometry 
of the lattice implies that the distribution is always centered and peaked around as indicated by the arrows. The 
additional peaks in the distribution reveal the discrete nature of the Si lattice. 

The respective distributions of exchange coupling between the same donor pairs in each ensemble is presented in 
Fig. |S1 where the arrows give the exchange value at the target relative position: J(Rt) ^ 0.29 meV. The results here 
are qualitatively different from the distance distributions in Fig. [S] since they are neither centered nor peaked at the 
target exchange value. Even for the smallest uncertainty radius of 1 nm in (a), the exchange distribution is peaked 
around J ~ 0, bearing no semblance to the inter-donor distance distributions. Increasing the uncertainty radius leads 
to a wider range of exchange values, with a more pronounced peak around the lowest J values. 

From the perspective of current QC fabrication efforts, ~ 1 nm accuracy in single P atom positioning has been 
recently demonstrated !34| . representing a major step towards the goal of obtaining a regular donor array embedded 
in single crystal Si. Distances and exchange coupling distributions consistent with such accuracy are presented in 
Figs. |3a) audita) respectively. The present calculations indicate that even such small deviations 1 nm) in the 
relative position of donor pairs with respect to perfectly aligned substitutional sites along [100] lead to order-of- 
magnitude changes in the exchange coupling, favoring J values. Severe limitations in controlling J would come 
from "hops" into different substitutional lattice sites. Therefore, precisely controlling of exchange gates in Si remains 
an open (and severe) challenge. 



C. Strained Si 



Uniaxial strain can be used to break the Si lattice symmetry and partially lift the degeneracy between the valleys, 
so that as few as two valleys make up the bottom of the conduction band. Then the sum over /i and v in Eq. H13|) 
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FIG. 6: Exchange distributions for the same relative position ensembles in Fig. |3 The arrow indicates the target situation. 
Contrary to the distance distributions, the exchange distributions are not centered or peaked around the target value. 



is much simplified, but the sinusoidal factor will still remain, so that care still has to be taken in controlling the 
donor exchange [39j. For example, if a uniaxial strain is applied along z direction, variation of relative donor position 
in the x-y plane does not lead to oscillation in the inter-donor exchange, as shown in Fig. [7| Strain is quantified 
here by a dimensionless parameter, x, defined in Ref. |89l |. On the other hand, atomic scale donor movement in a 
direction parallel to the strain direction z can still cause order-of-magnitude change in the donor exchange coupling, 
as is illustrated in Fig. |S| Therefore the exchange oscillation problem (leading to essentially random variations in the 
exchange coupling arising from dopant positioning within the Si unit cell) remains even in uniaxially strained Si since 
the system still has a two-valley degeneracy. These theoretical considerations ^] have recently been verified in a 
more sophisticated approximation which took into account higher-order Coulomb interaction corrections ^3 • 

Given the random variations in the exchange energy as a function of unintentional (and unavoidable at present) 
small 1 nm) variations in the positioning of P dopants in Si unit cell, as shown in Figs. I4I8I a careful rethinking 
of the exchange gate fabrication and control for Si QC architecture may be required. This is particularly true in 
view of the fundamental origin of the exchange oscillations, which arises from the quantum interference between 
the degenerate valleys in the silicon conduction band and therefore cannot be eliminated by simple fabrication tools 
(unless, of course, dopant positioning with sub-nanometer precision on the Si lattice somehow becomes available). 
One possibility is to eliminate the use of exchange gates and instead rely completely on hyperfine coupling and 
electron shuttling . Another possibility could be to go away completely from shallow donor bound states in Si and 
instead use SiGe quantum dot structures to confine single electrons, making the Si QC essentially identical to the 
proposed GaAs quantum dot QC |47| . Even in this case, however, the valley degeneracy problem may remain unless 
the quantum dot confinement potential is extremely weak (which would cause other problems such as the low lying 
orbital excited states potentially jeopardizing the two- level dynamics of the spin qubit). In any case, our results of the 
strong (essentially random) qubit-to-qubit variations in the exchange energy in the P-doped shallow-donor-based Si 
QC architecture compromises one of the main perceived advantages of Si QC (over, for example, the GaAs quantum 
dot QC), namely, the identical nature of each spin qubit (arising from all P shallow donor states in Si being identical 
by definition) does not seem to translate to similar exchange gate characteristics-in fact, the exchange gate turns out 
to be essentially random in nature (with wide variations in its strength peaking at zero). 

Comparing with Si, a question naturally arises about the GaAs quantum dot QC where the exchange would 
also exhibit random qubit-to-qubit variations arising from the (essentially trivial) random variations (±10%-20%) in 




FIG. 7: Calculated exchange coupling for in-plane displacements of the donors in the x-y plane for relaxed Si (x = 0, with 
diamond symbols) and strained Si (uniaxial strain along the z direction with x = ~20, with square symbols forming a straight 
line). The relative positions considered for the donor pairs consist of one donor at all possible lattice sites between two concentric 
circles of radii 90 A and 180 A with the other donor positioned at the center of the circles. The data points correspond to 
the exchange calculated at all relative positions considered. The solid line is J(R) for R along the [100] direction for x ~ —20. 




10 20 30 40 50 

Distancement from the reference site (Angstrom) 



FIG. 8; Donor electron exchange splitting in Si uniaxially strained along [001] direction (x = —20, corresponding to the strain 
in a Si quantum well grown in between relaxed Sio.8Geo.2 barriers). The two donors are approximately aligned along the [100] 
direction, with one of them displaced along the [010] direction (original figure appeared in Ref . l39tl . 
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quantum dot confinement sizes and inter-dot separations since it is impossible to fabricate identical quantum dots 
and/or to place them perfectly in an architecture. The quantum dot exchange variation problem (arising in the 
leading order from the e^^/" sensitivity of the exchange energy of the inter-dot separation R) is, however, less severe 
because the "random" exchange in this case will be peaked around a finite value in contrast to Si QC, where the 
most probable value of the inter-qubit exchange coupling is zero. We note, however, that in both cases one may 
have to carry out extensive characterization of the individual qubit-to-qubit coupling strength in order to operate the 
exchange gate, and such qubit characterization itself may turn out to be a hard control problem. We have recently 
proposed an experimental technique of spatially resolved micro-Raman spectroscopy as possible diagnostic tool to 
characterize local values of exchange coupling between individual spin qubits (i.e. the singlet-triplet energy splitting 
for individual pairs of induced shallow electron donor states) . 

IV. SUMMARY 

In summary, we have briefly reviewed two important issues related to donor-spin based quantum computing in 
silicon: quantum coherence of electron spins and exchange interaction among donor electrons. Our results show 
that the spin qubits based on donors in silicon have remarkable potential for very long quantum coherence times 
through isotopic purification. On the other hand, they also pose immense challenges in terms of precise nanostructure 
fabrications because of the degenerate nature of the silicon conduction band, which leads to essentially a random 
exchange gate coupling. Further studies of fabrication and innovative alternative approaches are imperative in order 
to fully realize the potential of this donor-based QC architecture. We have discussed an interesting dichotomy 
between donor-based Si and quantum dot-based GaAs QC architectures. In particular, isotopic purification, which 
can eliminate essentially all free nuclear spins (by eliminating ^^Si nuclei) in silicon (but not in GaAs), gives Si an 
enormous advantage of very long (T2 ^10— 100 ms) electron spin decoherence time compared with the corresponding 
GaAs quantum dot situation (T2 ~ 10 — 100 //s, three orders of magnitude shorter than in Si), making Si an ideal 
QC candidate material. However, the valley degeneracy of Si conduction band leads to demands for extremely precise 
dopant positioning (so that the inter-qubit exchange coupling is finite for most qubits) which will be difficult to 
achieve. It may be worthwhile in this context to consider alternative two-qubit gates, such as the dipolar gate or 
electron shuttling '4^ , for silicon quantum computation. 
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